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Abstract 

Background: A comparative analysis of the genomic and replication profiles of different geographical chikungunya 
virus (CHIKV) isolates of the East, Central and South African (ECSA) lineage was performed. 

Findings: Analysis of the data revealed the different growth kinetics for the different isolates. Deep genome 
sequencing analysis further revealed specific amino acid mutations in the viral nsP1, nsP3, nsP4, E1 and E2 proteins 
in the different isolates. Despite the difference in viral genomic profiles, the virus-induced ultrastructural changes 
within infected cells remained highly conserved among the different chikungunya virus isolates. 

Conclusions: These findings provide insights into the genomic and replication profiles of the re-emerging 
chikungunya virus isolates of the ECSA lineage. 
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Findings 

Chikungunya virus (CHIKV) was first isolated following 
an outbreak in Tanzania in the 1950s [1]. CHIKV infection 
typically induced an acute onset of painful syndrome as 
characterised by polyarthralgia, fever, asthenia, headache, 
myalgia and skin rashes [2,3]. 

CHIKV is an enveloped virus that possesses a single- 
stranded, positive-sense RNA genome. The 12Kb viral 
genome begins with an UTR at the 5'-terminal and 
followed by two-third of the coding regions for the non- 
structural proteins (nsPl-4), and the remaining one-third 
of the coding regions encodes for the structural proteins 
(Capsid-E3-E2-6K-El), and a ^-terminal poly-A-tail [4], 
The non-structural proteins are processed into four nsPs 
and responsible for different roles in virus replication. nsPl 
is important in viral RNA synthesis initiation and viral RNA 
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capping [5]. nsP2 possesses protease and RNA helicase ac- 
tivities [6,7]. nsP3 is composed of three domains and this 
viral protein is required for the formation and localization of 
replication complexes [7,8]. nsP4 possesses RNA-dependent 
RNA polymerase (RdRp) activity that is important for repli- 
cation and synthesis of the viral genome [8,9]. 

The urban transmission of the CHIKV involves the 
naive human host and the Aedes mosquito, with Aedes 
aegypti and Aedes albopictus as the main vectors [10]. 
Although the mortality rate is rather low, the epidemic 
potential of the virus is rather fearful due to the 
prolonged morbidity and the rapid spread of the virus. 
In 2004, the virus managed to catch fresh attention 
when it re-emerged and set its foot around the world. 
The first wave started in Kenya, the re-emergence of the 
virus was then reported in countries around the Indian 
Ocean [11,12], From 2006 onwards, more countries 
including India and Sri Lanka have also reported severe 
outbreaks of the CHIKV infections [13]. In addition, out- 
breaks of CHIKV were also documented in South East 



© 2013 Chen et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C6ntTcll Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Chen et al. Virology Journal 2013, 10:169 
http://www.virologyj.eom/content/1 0/1/1 69 



Page 2 of 8 



Asia, with countries such as Malaysia, Thailand and 
Singapore being affected [11]. 

Phylogenetic analysis of CHIKV genome revealed three 
lineages namely: Asian, West African (WA) and East, 
Central and South African (ECS A) [14]. Viral genome 
sequence analysis revealed that the CHIKV isolates from 
the recent outbreaks in the India Ocean and South East 
Asia were predominantly attributed by the ECSA lineage 
[12,15,16]. In addition, detailed phylogenetic analysis 
further identified two distinct sub-lineages in the ECSA 
clades between the outbreaks in Indian Ocean island 
and India [11]. The CHIKV outbreak in Indian Ocean, 
in particularly island of La Reunion, carries a distinct 
substitution of alanine for valine at the 226 position of 
the CHIKV El protein. This E1:A226V substitution has 
been suggested to allows a better adaptation of the virus 
in the Aedes albopictus as compared to classical vector, 
Aedes aegypti, which leads to a rapid dissemination of 
the virus in La Reunion [17]. In addition, the re-emergence 
CHIKV outbreak occurring in the South East Asia was 
determined to be more related to the Indian sub-lineage of 
ECSA as compared to the Indian Ocean [12]. 

Singapore had it first outbreak with CHIKV in the year 
of 2008 [18,19]. In 2008 alone, several outbreaks have 
occurred in Singapore with the first outbreak occurring 
in January 2008 and genome sequencing of the CHIKV 
isolate revealed that the virus does not harbour the 
distinct ELA226V mutation. The subsequent episodes of 
outbreaks in Singapore were caused by the CHIKV with 
the distinct ELA226V mutation [18]. Here, in this study, 
we are interested in the comparative analysis of the gen- 
omic sequences and replication profiles of the different 
Singapore CHIKV isolates and the La Reunion (IMT) 
CHIKV isolate that are belonging to ECSA lineage. The 
first CHIKV isolate, SGP007 is a wild type EL226A 
CHIKV isolated from patient in Singapore in 2008. The 
second CHIKV isolate SGP011, contains the ELA226V 
mutation and was also isolated from patient during the 
Singapore outbreak in 2008 [20]. The third isolate, La 
Reunion (IMT) was isolated from the outbreak in La 
Reunion that contains the ELA226V mutation [20]. All 
three isolates used were of low passage numbers upon 
isolation from patients. 

In the first part of the study, the virus infection was 
performed at a multiplicity of infection (MOI) of 10 and 
the replication kinetics of the three CHIKV isolates 
(SGP007, SGP011 and IMT) was analysed by negative- 
strand viral RNA production via quantitative RT-PCR 
(qRT-PCR) within Aedes albopictus (C6/36) cells and 
HeLa cells. Viral RNA was extracted from CHIKV- 
infected cell cultures using QIAamp viral RNA MINI kit 
(Qiagen, Hilden, Germany) and the qRT-PCR was 
performed using QuantiTect® Probe RT-PCR Kit 
(Qiagen, Hilden, Germany) with Applied Biosystems 



(ABI) 7900HT Fast Real-Time PCR System, using a 
TaqMan assay modified from Plaskon and colleagues 
[21]. In the C6/36 cells, we observed three distinct rate 
of growth kinetics for the three different CHIKV isolates. 
The ELA226V substitution in the El glycoprotein of the 
La Reunion isolate (IMT) allowed better adaptation in 
the C6/36 cells and resulted in a rapid increase in viral 
RNA copies as early as 6 hours post infection. The 
Singapore isolates (SGP007 and SGP011) started in- 
creasing in viral RNA copies from 12 hours post- 
infection. The SGP011 isolate, containing the ELA226V 
substitution in the El glycoprotein, showed higher cop- 
ies of viral RNA when compared to the SGP007 isolate, 
but still showed a poorer adaptability compared to the 
IMT isolate (Figure la). We then compared the replica- 
tion kinetics in HeLa cells, a susceptible human cell line. 
Generally, we saw the same trends as for the C6/36 cells, 
with a sharp increase in viral load in the IMT isolate 
when compared to the two Singapore isolates (Figure lb). 
Interestingly, it was noted that the infectious virus titres 
obtained by viral plaque assays further supported the no- 
tion that the IMT isolate produced a much higher infec- 
tious virus titre in both C6/36 and HeLa cells (Figure lc 
and d) when compared to both SGP007 and SGP011 
isolates. Notable, it was observed that SGP011 produced 
infectious viruses at an earlier time point of 12 hours 
post-infection when compared to SGP007 in C6/36 cells. 

We then proceeded to detect the presence of CHIKV 
antigens using flow cytometry in an effort to further 
define the virus replication kinetics. CHIKV-infected 
cells were fixed and permeabilised (BD Bioscience) be- 
fore being stained with an anti-alphavirus Ab (10 (ig/ml, 
Santa Cruz) for 30 min, followed by APC-conjugated 
goat anti-mouse IgG F(ab)2 Ab (10 (ig/ml, Invitrogen) 
for 30 min. Data was acquired in BD FACS Canto II 
(BD Biosciences) using BD FACSDiva software. In both 
C6/36 and HeLa cell lines, IMT and SGP011 isolates 
demonstrated similar infectivity kinetics and both iso- 
lates showed a significantly higher viral antigen staining 
compared to SGP007 (Figure le and f). This suggested 
that both IMT and SGP011 exhibited higher virus infec- 
tivity and viral antigen production in the infected cells 
when compared to SGP007 isolate. 

To better understand if the differences in replication 
kinetics of the different CHIKV isolates were associated 
with specific mutation(s) in the viral genomes, viral 
RNA was isolated from each CHIKV isolate and reverse- 
transcribed into cDNA and sent for sequencing. The 
three cDNA libraries were sequenced as a multiplex in a 
single lane using the Illumina HiSeq2000 (Next Gene- 
ration Sequencing Core facility, Genomic Institute of 
Singapore) at a coverage of 67436 x, representing about 
20 million reads per sample. De novo assembly and map- 
ping of each sequence was performed using the velvet 
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Figure 1 Replication kinetics of CHIKV isolates (IMTSSA, SGP011 and SGP007) in infected C6/36 and HeLa cells, a and b) Negative 



strands viral RNA quantification with qRT-PCR of CHIKV infected C6/36 and HeLa cells over a period of five days post-infection, c and d) Infectious 
virion formation from CHIKV infected C6/36 and HeLa cells were analyzed by plaque formation assay over a period of five days post-infection and 
expressed as log PFU/ml. e and f) FACS quantification of CHIKV antigen in infected C6/36 and HeLa cells over a period of five days post-infection. 
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assembler and bowtie enabled by Pipeline Pilot, with hash 
size = 31 and insert size = 150. All three sequences 
seemed to be missing 20-50 bp from the 3' and 5' ends 
when compared to the prototype isolate (S27 African 
prototype). In order to perform end-repair, 100 000 reads 
were mapped with bowtie to the prototype isolate, and 
some of the end sequences were recovered. The libraries 
were compared with the ERCC spike-in RNA controls 
(Invitrogen), and all three libraries had a > 0.99 correlation 
between normalised read counts and the control concen- 
trations in a log scale. 

When compared to the CHIKV S27 African prototype 
isolate [22], we identified seven nucleotide (nt) substitu- 
tions within the genome of the SGP007 CHIKV isolate 
that resulted in amino acid (aa) changes (Additional file 1: 
Table SI). S27 African prototype and SGP007 CHIKV 
isolates had identical amino acid residues at various posi- 
tions (nsP3 aa524, nsP4 aa75 and El aa226). The El 
A226V mutation suggested a closer identity between these 
two isolates. However, SGP007 CHIKV obtained an extra 
unique mutation at various positions (nsPl aa60, nsP3 
aal37, nsP4 aa563 and E2 aal78). It remains unknown 
whether these mutations will support the unique feature 
of SGP007 CHIKV in terms of virus replication. We did 
not identify any unique mutations within the nsP2, Capsid, 
E3 or 6K protein sequences. 

Till now, there is limited information available regard- 
ing the relationship between the unique amino acid resi- 
dues and the functions of alphavirus nsPs. Previous 
studies have suggested the role of El A226V mutation in 
the fitness of CHIKV in Aedes aegypti and Aedes 
albopictus mosquitoes [23]. Further studies have also 
suggested the role of E2 L210Q mutation in the sequen- 
tial adaptive mutations in CHIKV and enhancement of 
efficient vector switching in nature [24]. More amino 
acid positions within the CHIKV structural proteins have 
been demonstrated to show positive selection effect in 
virus circulation and persistence in endemic areas [25]. 
Here, we would like to postulate the role and effect of 
amino acid mutations in the nsPs on virus replication 
efficiency. 

We first identified one unique mutation in SGP007 
(nsP3 A137V) located within the macro domain region 
(Additional file 1: Table SI). The first 160 amino acids of 
the N-terminal region of nsP3 are known as the macro 
domain, which possesses ADP-ribose l"-phosphate 
phosphatase activity. In addition, a study using Sindbis 
virus suggested an important role of the nsP3 macro 
domain in virus replication [26]. However, structural 
mapping suggested that this amino acid is not positioned 
within the enzymatic groove of the protein (data not 
shown). Therefore, it is unlikely that mutations at this 
position would alter the phosphatase activity and further 
studies would be needed to investigate the importance 



of this mutation (nsP3 A137V) in CHIKV virus 
replication. 

We further suggested a possible role for the opal stop 
codon that located at the C-terminal of the nsP3 
proteins in enhancing CHIKV virus replication in vitro. 
Several alphaviruses express two different isoforms of 
nsP3, depending on the presence of an opal stop codon or 
an open reading frame upstream of nsP4 (7 amino acids 
before the C-terminal of the nsP3) [27,28]. Expression of 
CHIKV nsP3 in different isoforms with different proper- 
ties or functions remains to be further characterized. 
CHIKV S27 African prototype and SGP007 have an iden- 
tical amino acid residue at nsP3 position 524 (arginine), 
but in IMT and SGP011, this position encodes a stop 
codon instead. In vitro infection results showed that 
SGP007 resulted in lower virus replication efficiency when 
compared to IMT and SGP011 (Figure 1). It was shown 
previously that the Opal stop codon at this particular 
position improved O'Nyong-Nyong virus infectivity, rela- 
tive to viruses encoding an arginine at the corresponding 
position [29]. In addition, a study using Semliki Forest 
virus also showed improved nsP3 protein stability after 
removing the nsP3 C-terminal region [28]. Therefore, our 
results are consistent with previous findings from other 
alphaviruses that CHIKV isolates possessing an opal stop 
codon upstream of nsP4 could improve nsP3 protein sta- 
bility and virus infectivity. This in turn could enhance 
virus replication. In addition, the T75A amino acid residue 
at the nsP4 N-terminal region was found to correlate with 
the nsP3 C-terminal region in regulating CHIKV viral 
replication. 

One last unique mutation identified from CHIKV 
SGP007 was the I563T amino acid residue at nsP4 
located within the RdRp region. Within the structural 
region, one unique mutation identified from CHIKV 
SGP007 was the R178H amino acid residue at E2 glyco- 
protein and is located within the E2 domain B region. 
However, the precise role of the differences in these 
amino acids in the various CHIKV isolates in virus repli- 
cation remains undefined. 

To further address the issues if the specific amino 
acids differences between the different CHIKV isolates 
will affect the formation and morphogenesis of virus - 
induced structures during CHIKV replication within 
cells, we utilised transmission electron microscope 
(TEM) as a tool to decipher these replication processes 
within the infected cells. The replication process of 
CHIKV has not been well-defined and many studies are 
still based on the assumption on other established 
alphaviruses. In this part of the study, the C6/36 cells 
were infected with CHIKV at the MOI of 10 and mor- 
phological changes were observed and documented 
under TEM at different time points post-infection. Simi- 
larly, the HeLa cells were infected with CHIKV at the 
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same MOI and processed for the ultrastructural analysis 
at different time points post-infection with the different 
isolates of CHIKV. 

In the mocked-infected cells, normal morphology of 
the cells was observed and recorded (Figure 2a). The 
C6/36 cells were characterized by large nuclei, occupying 
two-thirds of the cells and the other cellular organelles 
were also observed. At higher magnification, stacks of 
Golgi apparatus, smooth and rough endoplasmic reticu- 
lar and mitochondria were clearly identified and docu- 
mented (Figure 2a). 

Generally, upon infection with the different CHIKV iso- 
lates, the virus induced structures and morphogenesis 
observed were highly conserved. The first distinct viral- 
induced structures observed, were the presence of the cy- 
topathic vacuoles type I (CPV I) complexes. This occurred 
as early as 6 to 8 hours post-infection. These distinct vacu- 
oles measured to be approximately 800 nm to 2000 nm in 
diameter and are characterized by the presence of mem- 
branous spheres lining on the vacuoles (Figure 2b). These 
membranous spheres (Figure 2b, arrows) measured to be 
60 nm to 80 nm in sizes and were seen to be neatly 
arranged on the vacuoles. These vacuoles were also docu- 
mented in cells infected with other alphaviruses and these 
vacuoles were proven to be the replication complexes 
where viral RNA synthesis occur and were seen in the 
early infection [30]. Other studies also showed the pres- 
ence of the viral replicases in CPV I complexes of other 
alphaviruses. The staining of the cellular markers for the 
cytoplasmic endosome and lysosome suggested the origin 
of these CPV I complexes [31,32]. As the infection 
progressed to 10 hours post-infection, the degenerated 
forms of CPV I can be observed. The membranous 
spheres within the CPV I complexes are no longer visible 
and are replaced by vacuoles containing bags of membra- 
nous whorls (data not shown). 

The appearance of the CPV I complex at 8 hours post- 
infection was accompanied by the emergence of the 
mature virus particles budding at the plasma membrane 
(Figure 2c). The viral nucleocapsids of the CHIKV mea- 
sured approximately 30 nm to 40 nm in size can be seen 
to be lining in cytoplasm, underneath the plasma mem- 
brane (Figure 2ci). The viral nucleocapsids were seen to be 
expelled out of the cytoplasm through blebbing and 
pinching off the plasma membrane (Figure 2d), acquiring 
a portion of the membrane as its envelope (Figure 2di). 
The full mature virus was measured to be approximately 
60 nm to 70 nm in size. The progression of the infection 
then saw an array of neatly arranged mature virus lining in 
rows on the peripheral of cells and in between the 
junction of two adjacent cells. This mode of virus release 
from the infected cells via budding through the plasma 
membrane occurred at early stages of infection and 
marked the first mode of viral release. 



As the infection progressed, the dominant viral- 
induced structure was replaced by the cytopathic vacuoles 
type II (CPV II) complexes. This phenomenon occurred as 
early as 24 hours post-infection. Unlike the CPV I, these 
CPV II complexes were characterized by smaller vacuoles, 
measuring 300 nm to 400 nm in size (Figure 2e). These 
complexes were often observed in large number and clus- 
ter near the plasma membrane (Figure 2f). The cyto- 
plasmic side of the CPV II vacuoles was surrounded by 
neatly arranged electron-dense viral nucleocapsids meas- 
uring 30 nm to 40 nm in size. These nucleocapsids can be 
observed to be budding into the CPV II complexes, and 
were thought to mature in these CPV II vacuoles. Subse- 
quently, these CPV II vacuoles can be seen to contain full 
matured enveloped CHIKV viruses (Figure 2f, arrow), and 
the fusion of these vacuoles with the plasma membrane to 
release the matured virus out of the cell (Figure 2g). As 
such, it can be postulated that the CPV II may function as 
exocytosis vesicles or vacuoles for the release of the 
CHIKV particles at late infection. In addition, crystalline 
array of virus particles accumulating on the surface of the 
infected cells can be seen (Figure 2h). This mode of virus 
release through exocytosis of the CPV II complex with the 
fusion with the plasma membrane marked the second 
mode of viral release during the late stages of infection 
and may provide a means of rapid release of the mature 
viruses to the exterior of cell before cell death sets in. 

CHIKV infection was also tracked in the HeLa cells 
with CHIKV infection. The mock-infected HeLa were 
characterized with typical cellular organelles as shown 
in Figure 3a. Generally, for the CHIKV-infected cells, 
the appearance of the distinct viral-induced structures 
remained similar to that of the infected C6/36 cells. 
Interestingly, CPV I complexes were no longer the pre- 
dominant viral-induced structures observed at 6 hours 
post-infection, instead extensive membranous spheres 
were also observed to be lining on the plasma mem- 
brane (Figure 3b, boxed). These membranous spheres 
were measured to be approximately 60 nm to 80 nm 
and were structurally identical to the membranous 
spheres in the CPV I complexes (Figure 2b, arrows). It 
was presumed that these membranous spheres were the 
RNA replication site for the viral synthesis. The forma- 
tion of these replication complexes was also noted to be 
the initial replication complexes in the vertebrate cells 
but not the invertebrates cells in other alphaviruses 
[33]. As infection progressed, matured viruses were seen 
to be budding off the plasma membrane as early as 8 
hours post-infection (data not shown). The formation of 
extensive CPV II complexes can be noted as early as 12 
hours post-infection. Rows of viral nucleocapsids can 
also be seen to be budding into the vacuoles extensively 
at the different planes of the cut section (Figure 3c). 
Similar to the infection in the C6/36 cells, the CPV II 
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(See figure on previous page.) 

Figure 2 Replication of CHIKV (IMT isolate) in C6/36 cells, (a) Ultrastructural analysis of mock-infected C6/36 cells. At higher magnification, 
Golgi apparatus (GA), endoplasmic reticulum (ER) and mitochondria (M), can be observed. The bar corresponds to 0.5 urn. (b) Presence of 
Cytopathic Vacuoles (CPV) I in early infection (8 h post-infection). A typical CPV I complex (appx 800 nm) can be observed. CPV I is characterized 
by membranous spherules (arrows) regularly arranged from the interior surface of the complex. The bar corresponds to 0.2 urn. (c) Virus release 
by budding at the plasma membrane during early infection. CHIKV particle are observed to be budding from the plasma membrane (Pm). The 
bar corresponds to 0.2 urn. (c) i Enlarged of (c) showing the lining of the nucleocapsids (arrowheads) underneath the plasma membrane (Pm) 
of the cell. Matured virus (arrow) is observed to be budding from the peripheral of cells membrane, (d) Budding of CHIKV occurred extensively 
at the plasma membrane of the cells. The bar corresponds to 0.2 urn. (d) i Enlarged of (d) monitored the different stages of budding of CHIKV 
particles from the plasma membrane (arrows), (e) Ultrastructural analysis of late infection with CHIKV in C6/36 cells (Day 1 post infection). CPV II 
formations are observed. The viral nucleocapsids (arrowheads) are present at the cytoplasmic surface of CPV II complexes. The bar corresponds 
to 0.2 urn. (f) Several CPV II complexes are found in cluster near the plasma membrane (Pm) with the noticed of mature virus (arrows) in the 
complex, (g) and (h) At the late stage, CHIKV can be seen to be released from cells via exocytosis (arrows). Both the bar of g and h correspond 
to 0.2 urn. 



complexes were often present in cluster and located in 
close association to the plasma membrane. However, the 
number of CPV II complexes formed in HeLa cells was 
much more extensive than the infection in the C6/36 cells 
(Figure 3c). The rapid maturation and proliferation of the 
CHIKV virus might lead to the failure of proper exocytosis 
of the matured virus through fusion of the plasma 
membrane and resulted in the massive accumulation of 
matured viruses within the infected HeLa cells (Figure 3d). 



This may ultimately lead to the cell death of the CHIKV- 
infected Hela cells. However, it is also important to note 
that Krejbich-Trotot and workers (2011) has recently illus- 
trated that CHIKV could also mobilize the apoptotic 
machinery by extrinsic and intrinsic pathways in Hela 
cells, a mechanism that could also be used to escape 
macrophages and viral host defences in vivo [34]. 

The re-emergence of CHIKV since 2004 has brought 
in much attention to this once neglected virus. This 




Figure 3 Replication of CHIKV (IMT isolate) in HeLa cells, (a) Ultrastructural analysis of mock-infected HeLa cells. Under high magnification, 
the cellular organelles such as the Golgi apparatus (GA) and mitochondria (M), can be seen. The bar corresponds to 0.2 urn. (b) Infection of CHIKV 
in HeLa cells. Viral-induced membranous spheres are observed to associate with the plasma membrane [boxed]. The bar corresponds to 0.2 urn. 
(c) CPV II complexes were detected at 12 hours post infection. The viral nucleocapsids (arrows) were observed to be budding in the CPV II 
complexes. The bar of c corresponds to 0.5 urn. (d) Massive accumulation of mature CHIKV (arrows) in infected cells at late infection. The bar 
corresponds to 0.5 urn. 
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study has documented the differences in the replication 
kinetics of the different CHIKV isolates within the ECSA 
lineage. Despite the differences in the replication kinetics, 
the virus induced ultrastructural changes and virus 
morphogenesis remains highly conserved between the 
different CHIKV isolates of the ECSA lineage. 

Additional file 



Additional file 1: Table SI. Amino acid differences between the CHIKV 
S27 sequence (template) and three other CHIKV isolates (IMT, SGP007 
and SGP011). 
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